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Abstract. In this work we consider the problem of the approximate hedging of a contingent 
claim in minimum mean square deviation criterion. A theorem on martingale representation in 
the case of discrete time and an application of obtained result for semi-continous market model 
are given. 

Keywords: Hedging, contingent claim, risk neutral martingale measure, martingale represen- 
tation. 


1. Introduction 


The activity of a stock market takes place usually in discrete time. Unfortunately such markets 
with discrete time are in general incomplete and so super-hedging a contingent claim requires usually 
an initial price two great, which is not acceptable in practice. 

The purpose of this work is to propose a simple method for approximate hedging a contingent 
claim or an option in minimum mean square deviation criterion. 


Financial market model with discrete time: 


Without loss of generality let us consider a market model described by a sequence of random 
WECLOIS® {8a 7 Opel ayey fd R?, which are discounted stock prices defined on the same 
probability space {Q, S, P} with {F,, n = 0,1,...,.N} being a sequence of increasing sigma- 
algebras of information available up to the time n, whereas “risk free ” asset chosen as a numeraire 
so 1. 

A Fyx-measurable random variable H is called a contingent claim (in the case of a standard call 
option H = max(S, — K,0), K is strike price. 
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Trading strategy: 


A sequence of random vectors of d-dimension y = (yn, 2 = 1,2,..., N) with yn, = (71, -72,..., 
42)? (A? denotes the transpose of matrix A ), where ae is the number of securities of type 7 kept by 
the investor in the interval [n — 1, n) and 4, is F,-1 -measurable (based on the information available 
up to the time n — 1), then {7,} is said to be predictable and is called portfolio or trading strategy . 


Assumptions: 


Suppose that the following conditions are satisfied: 
i) AS, = Sp — Sn-1, H € Lo(P), n=0,1,...,N. 
ii) Trading strategy + is self-financing, i.e. S7_,yn-1 = S7_,%n or equivalently S7_, Ayn = 0 
for allt 1/2) :.9n8V¥- 
Intuitively, this means that the portfolio is always rearranged in such a way its present value 
is preserved. 
iii) The market is of free arbitrage, that means there is no trading strategy such that 7 So := 
71-50 <0, yw.Sn > 0, Pyn.Sy > 0} > 0. 
This means that with such trading strategy one need not an initial capital, but can get some profit and 
this occurs usually as the asset {Sn} is not rationally priced. 


Let us consider 
d 


N 
Gu (1) = >> e-ASk with AS, = 5° As}. 
k=1 j=) 
This quantity is called the gain of the strategy y . 
The problem is to find a constant c and y = (yn, n = 1,2,..., N) so that 


Ep(H —c- Gy(7))? > min. (1) 


Problem (1) have been investigated by several authors such as H.folmer, M.Schweiser, M.Schal, 
M.L.Nechaev with d = 1. However, the solution of problem (1) is very complicated and difficult for 
application if {S,} is not a {F,,}-martingale under P, even for d = 1. 

By the fundamental theorem of financial mathematics, since the market is of free arbitrage, there 
exists a probability measure Q ~ P such that under Q {S,} is an {F,,}-martingale, i.e. Eg(Sn| Pn) = 
Sn-1 and the measure Q is called risk neutral martingale probability measure . 

We try to find c and 7¥ so that 


Eo(H —c— Gn(7))* > min over +. (2) 


Definition 1. (7%) = (y%(c)) minimizing the expectation in (1.2) is called Q- optimal strategy in the 
minimum mean square deviation (MMSD) criterion corresponding to the initial capital c. 

The solution of this problem is very simple and the construction of the Q-optimal strategy is 
easy to implement in practice. 

Notice that if Ly = dQ/dP then 


EQ(H — c— Gy(7))? = Ep[((H — c- Gw)?Ly] 


can be considered as an weighted expectation under P of (H — c— Gy)? with the weight Ly. This 
is similar to the pricing asset based on a risk neutral martingale measure Q. 
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In this work we give a solution of the problem (2) and a theorem on martingale representation 
in the case of discrete time. 

It is worth to notice that the authors M.Schweiser, M.Schal, M.L.Nechaev considered only the 
problem (1) with S, of one-dimension and M.Schweiser need the additional assumptions that {S,} 
satisfies non-degeneracy condition in the sense that there exists a constant 6 in (0,1) such that 


(E[ASp|Fr-1})? < 6E[(ASp)?|Fn-1] — P-a.s. for all n = 1,2,..., N. 
and the trading strategies +,,’s satisfy : 
ElynAS,]? < 0, 


while in this article {S,,} is of d-dimension and we need not the preceding assumptions. 

The organization of this article is as follows: 

The solution of the problem (2) is fulfilled in paragraph 2.(Theorem 1) and a theorem on the 
representation of a martingale in terms of the differences AS, (Theorem 2) will be also given (the 
representation is similar to the one of a martingale adapted to a Wiener filter in the case of continuous 
time). 

Some examples are given in paragraph 3. 

The semi-continuous model, a type of discretization of diffusion model, is investigated in para- 
graph 4. 


2. Finding the optimal portfolio 


Notation. Let Q be a probability measure such that Q is equivalent to P and under Q {S,, n = 
1,2,..., NM} is an integrable square martingale and let us denote En(X) = Eg(X|Fn), Hn = 
H, H, = Eo(Al|Fn) = En(H); Vatn-1(X) = [Covn-1(Xi, X;)] denotes the conditional variance 
matrix of random vector X when F;,-, is given, I’ is the family of all predictable strategies +. 


Theorem 1. /f {S,} is an {F,,}-martingale under Q then 


Eq(H — Ho — Gw(7*))? = min{ E9(H — c— Gw())?: €T}, (3) 
where ‘y~ is a solution of the following equation system: 
[Varg=1(AS,) 425 = Enay((AH,AS,) * -P- as. (4) 


Proof. At first let us notice that the right side of (3) is finite. In fact, with y, = 1 for all n, we have 


2 
N d 
Eo(H -c-Gy(7))?=Eg|H-c->_>_ wes | < 00. 
n=1 j=1 ) 
Furthermore, we shall prove that y*AS,, is integrable square under Q. 

Recall that (see [Appendix A]) if Y, X, Xo,..., Xz are d+1 integrable square random variables 
with E(Y) = E(X;) =--- = E(Xqg) = 0 and if Y = b)X1+boXo+--+ +bgXq is the optimal linear 
predictor of Y on the basis of X,, Xo,...,.Xq then the vector b = (6, b2,...,6g)7 is the solution of 
the following equations system : 

Var(X)b = E(YX), (5) 
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and as Var(X) is non-degenerated b is defined by 


b = [Var(X)] 'E(YX), (6) 
and in all cases 
bTE(YX) < E(Y?), (7) 
where X = (X1, Xo,..., Xe). . 
Furthermore, 
Y -Y1LX;j, ie. E[X(¥ —Y)] =0, i=1,...,k. (8) 


Applying the above results to the problem of conditional linear prediction of AH, on the basis 
of AS}, AS?,..., AS? as F, is given we obtain from (5) the formula (4) defining the regression 
coefficient vector y*. On the other hand we have from (5) and (7): 


E( 7,7 ASn)? = EEn-i(7,) ASp ASE 9g") = E(yg7 Vatn—1(ASn) 90) 
= E(¥2En-1(AHnASn)) < E(AHn)? < 00. 


With the above remarks we can consider only, with no loss of generality, trading strategies +, such 
that 


En-1(%nASn)* < 00. 
We have: 
Hn = Ho + A, +::: + AHn 
and 


En-1(AHn — Y2ASn)? = En-1(ANn)? — 297 En-1((AHn ASn) + YF En-1(ASn AS? )Yn- 


This expression takes the minimum value when 7, = 7;.. 
Furthermore, since {H, — c — G,(7)} is an {F,,}- integrable square martingale under Q, 


N 2 
Eq(Hw — c— Gn(7))? = Eq | Ho - ¢- (An — inASn) 
n=] 
"oN 2 
= (Hp —c)? + Eg (AH, —- ic) 
n=] 
N 
= (Ho —c)* + > Eo(AHn — %nASn)? (for AHn — YynASn being a martingale difference) 
n=1 


N 

= (Ho — c)? + Eg ) > En-1(AHn — YnASn)? 
n=) 
N 

> (Ho -— c)? + Eq > En-1(AHn - 7345p)" 


n=1 
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N 
= (Ho - c)* + Eg )-(AHn - ¥,ASn)? 
n=] 


N 2 


> (AHn — ¥ASn) 


n=) 


= (Ho -¢)? + Eg 


> Eg(Hn -— Ho - Gn(7"))?. 


So Eg(Hn — c- Gy(7))? > Eo(Hn — Ho — Gn(y"))? and the inequality becomes the equality if 
o> Hoard y= 4: 


3. Martingale representation theorem 


Theorem 2. Let {H,, n= 0,1,2,...}, {Sp, n =0,1,2,...} be arbitrary integrable square random 
variables defined on the same probability space {Q,3,P}, FS = o(So,...,Sn). Denoting by 
M1(S, P) the set of probability measures Q such that Q ~ P and that {S,,} is {F2} integrable square 
martingale under Q, then if F = \/°, FS, Hn, Sn € Lo(Q) and if {H,} is also a martingale under 
Q we have: 


1. Hn = Ho t+ >) AS + Cr a.s., (9) 
k=] 
where {Cn} is a {>} —Q-martingale orthogonal to the martingale {Sn}, ie. En—\((AC,AS_) = 0, 
for all n =0,1,2,.., whereas {yn} is {F3_,}- predictable. 


2. Hn = Ho+ >> yf ASk = Hot+Ga(y) Pas. (10) 
k=1 
for all n finite iff the set T1(S, P) consists of only one element. 
Proof. According to the proof of Theorem 1, Putting 


ACy = AHk — ,7ASk, Cn = D> ACK, Co = 0, (11) 
k=l 
then AC, LAS,, by (8). 
Taking summation of (11) we obtain (9). 
The conclusion 2 follows from the fundamental theorem of financial mathematics. 


Remark 3.1. By the fundamental theorem of financial mathematics a security market has no arbitrage 
opportunity and is complete iff II(S,P) consists of the only element and in this case we have (10) 
with yy defined by (4). Furthermore, in this case the conditional probability distribution of S,, given 
Fe | concentrates at most d + 1 points of R? (see [2], [3]), in particular for d = 1, with exception of 
binomial or generalized binomial market models (see [2}, [7]), other models are incomplete. 


Remark 3.2. We can choose the risk neutral martingale probability measure Q so that Q has minimum 
entropy in II(.S, P) as in [2] or Q is near P as much as possible. 
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Example 1. Let us consider a stock with the discounted price Sp at ¢ = 0, S, at t = 1, where 
4S9/3 with prob. p), 
Si = 4 So with prob. pe, p1,p2,p3>0, pi +pe2+p3=1 
5509/6 with prob. p3. 
Suppose that there is an option on the above stock with the maturity at £ = 1 and with strike price 
K = So. We shall show that there are several probability measures Q ~ P such that {So, Si} is, 
under Q, a martingale or equivalently Eg(AS,) = 0. 

In fact, suppose that Q is a probability measure such that under Q S takes the values 
4So/3, So, 2509/3 with positive probability qi, ¢2, q3 respectively. Then Eg(AS,;) = 0 © 
So(q1/3 — 93/6) = 0 <> 2g; = q3, so Q is defined by (q1, 1 — 3q1, 291), O< qm < 1/3. 

In the above market, the payoff of the option is 

H = (Si — K)+4 = (AS))4 = max(AS;, 0). 
It is easy to get an Q-optimal portfolio 
7° = EQ[HASi)/Eq(451)? = 2/3, EQ(H) = 4150/3, 
Eg|H — Eq(H) - 7°AS;]? = 9 55(1 — 3q1)/9 — 0 as qi > 1/3. 
However we can not choose q, = 1/3, because g = (1/3, 0, 2/3) is not equivalent to P. It is better 


to choose g1 = 1/3 and 0 < q; < 1/3. 
Example 2. Let us consider a market with one risky asset defined by : 


n 
So = So[[ Zi: or ron = Sriens n= ED. mags 

i=] 
where Z}, Z2,..., Zn are the sequence of i.i.d. random variables taking the values in the set {2 = 
{d,,dg,...,dyz) and P(Z; = d,) =p, > 0, K=1,2,..., M. It is obvious that a probability measure 
Q is equivalent to P and under Q {S,} is a martingale if and only if Q{Z; = dx) = qm > 0, k = 
1,2,...,M and Eg(Z;) =1, ie. 

qd) + godo2+-++:+qudy =1. 
Let us recall the integral Hellinger of two measure Q and P defined on some measurable space 
{N°, F}: 

H(P,Q)= [ (ap.aqy'”? 
ge 

In our case we have 


H(P,Q) = > {P(Z, = da, Ze = dia,..., Zv = din)"Q(Z, = di, Z2 = dia,..., Zn = din)'? 
= So {pigin pi2gie-.-pingin}'? 


where the summation is extended over all dj, di2,..., din inQ or over all 71, i2,...,i2n in {1,2,..., M}. 
Therefore 


(4 ‘ 
H(P,Q) = Sinai | ; 


. i=] 


We can defime a distance between P and Q by 


IQ — Pil? = 2(1 - H(P,Q)). 
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Then we want to choose Q* in II(S, P) so that ||Q* — P|| = inf{||Q — P|: Q € II(S, P)} by solving 
the following programming problem: 


ual 2 

lee } 
yo 77g, — max 
i=] 


with the constraints : 
i) gid) + qoad2+ ---+qQmdm =1. 
li) a +go+--> +a =1. 
iii) Qi, 92,--+. Oa > 0. 
Giving pj, p2,.--, Day we can obtain a numerical solution of the above programming problem. It is 


possible that the above problem has not a solution. However, we can replace the condition (3) by the 
condition 


iii’) Qi, G2;+++3 Wd 2 0, 
then the problem has always the solution g* = (q7, ¢3,.--, 94) and we can choose the probabilities 
Qi; G2)--+,qQn¢ > O are sufficiently near to qj, q9,---: 944 - 


4. Semi-continuous market model (discrete in time continuous in state) 


Let us consider a financial market with two assets: 
+ Free risk asset {B,, n = 0,1,..., N} with dynamics 


n \ 
By i= exp (>: rs} Sis sek. (12) 


k=1 
+ Risky asset {S,, n = 0,1,..., N} with dynamics 


Sn = So exp pace + o(S-s)al (13) 
k=1 
where {gn, 2 =0,1,..., N} is a sequence of i.i.d. normal random variable N(0, 1). It follows from 
(13) that 
Sp = Sp-1 exp((Sp-1) + 0(S, -1)9n); (14) 


where Sp is given and 44(S,-1) := a(Sp-1) — o7(Sp—1)/2, with a(x), o(z) being some functions 
defined on [0, 00) . ; 
The discounted price of risky asset S, = S,,/B, is equal to 


Sn = Soexp pazcey —Tre+ o(Se-s)a , (15) 


k=1 
We try to find a martingale measure Q for this model. 
It is easy to see that Ep(exp(Agx)) = exp(A?/2), for gz ~ N(0, 1), hence 


Eexp (Stasi - px Sia%/l =1 (16) 


k=1 
for all random variable @;(S%-1) . 
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Therefore, putting 


Ln = exp (Stasi) = pxsi-1)/2) n= Leuk (17) 
k=1 
and if Q is a measure such that dQ = LxdP then Q is also a probability measure. Furthermore, 
= = exp(p(Sn-1) = Faget o(Sn-1)9n)- (18) 
Denoting by £E°, E expectation operations corresponding to PP, QO, 
E,(.) = E[(.)|F5] and choosing 


then it is easy to see that 
En-1(Sn/Sn-1| = E°[LnSn/Sn-1|Fo]/Ln-1 = 1 
which implies that {S,} is a martingale under Q. 
Furthermore, under Q, S;,, can be represented in the form 
Sn = Sn-1 exp((u"(Sn-1) + 0(Sn-1) 95). (20) 


Where :*(Sp-1) = Tn — 07 (Sn-1)/2, 9% = —Bn + 9n is Gaussian N(0, 1). It is not easy to show the 
structure of II(S, P) for this model. 

We can choose a such probability measure & or the weight function Ly to find a Q- optimal 
portfolio. 
Remark 4.3. The model (12), (13) is a type of discretization of the following diffusion model: 

Let us consider a financial market with continuous time consisting of two assets: 


+Free risk asset: 
t 
By = exp ¢ r(u)du) : (21) 
fn) 


+Risky asset: dS; = S,[a(S:)dt + of(S;:)dWt], So is given, where 
a(.), o(.) : (0,00) — R such that za(x), zo(x) are Lipschitz. It is obvious that 


S; = exp { fais.) — o*(Sy)/2}du + i: a(s.)aw.} » Wiehe Te. (22) 
0 4] 
Putting 
p(S) = a(S) - o°(S)/2, (23) 


and dividing [0,7] into N intervals by the equidistant dividing points 0,A, 2A,..., NA with 
N =T/A sufficiently great, it follows from (21), (22) that 


nd nd | 
p(S,,)du + | o(Sy)dW,, 

meh. (n—-1)A 
= Sin-1)a exP{H(S(n-1ya)A + (Stn—1ya) [Waa — Wn-1al} 
© Sin-1)a exP{u(Sin—1)a)O + o(S(n—1)a) A? 9n} 


Sna = S(n-1)a €XP 
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with gn = [Wana — Wen-1yal/Al/?, n = 1,..., N, being a sequence of the i.i.d. normal random 
variables of the law N(0, 1), so we obtain the model : 

Sra = Sina XP{H(Shr_1ya)d + (Shr—-1ya) 4/7 9n}- (24) 
Similarly we have 

Bra = Bian-1ya &xP(T(n-1)a4).- (25) 

According to (21), the discounted price of the stock 5; is 

= S ¢ ) 

Ss; = B, = So exp t[ (u(Su) — ruldy + i o(Su)aW¥u p : (26) 

0 0 


By Theorem Girsanov, the unique probability measure Q under which {5,, F5, Q} is a martingale 
is defined by 
s = wie 
(aQ/aP Fe = exp (faa, - 5 | Bed, := Lew), (27) 
) 


where 
3, = ((a(Ss) = 75) 
o_ — gxre ys,  4t 
o(S;) 
and (dQ/dP)|F? denotes the Radon-Nikodym derivative of Q w.r.t. P limited on F, Furthermore, 
under Q 


t 


0 
is a Wiener process. It is obvious that LT can be approximated by 


N ‘ 
Ly := exp (>> A, A'!?q, — Aot/2 (28) 
k=l 
where (a(S 
_— Q\O(n-1)A) — Tnd 929 
Zs a(S(n-1)a) ae) 


Therefore the weight function (25) is approximate to Radon-Nikodym derivative of the risk unique 
neutral martingale measure Q w.rt. P and Q is used to price derivatives of the market. 


Remark 4.4. In the market model Black- Scholes we have Ly = Lr. We want to show now that for 
the weight function (28) 

Eg(H — Ho - Gn(y"))? +0 as N 4 co or A 0. 
where -y* is Q-optimal trading strategy. 


Proposition. Suppose that H = H(Sr) is a integrable square discounted contingent claim. Then 

Eo(H — Ho — Gn(7*))? +0 as N - cco or A 0, (30) 
provided a, r and o are constant (in this case the model (21), (22) is the model Black-Scholes ). 
Prcof’. It is well known (see[4], [5]) that for the model of complete market (21), (22) there exists a trad- 
ing strategy » = (ym =~ gf(t,S(t)), O = %t = = T), hedging 
H, where vy : [0, T] x (0,00) — RF is continuously derivable in ¢ and S, such that 


T _— 
H(Sr) = Ho + [ yidS(t) as. 
0 
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On the other hand we have 


N 2 
EQy (w - Ho- dU va-yaddas] 
\ k=1 


N 2 
Son (« ee > yuna 
k=1 
P 2 


T ‘ N 7 \ 
= Eq ( yrds (t) — 2, (n-aASin-a } Ly/Lr 
k=1 


T N 2 

= Eq (| yrd S(t) — Y#u-naASn-a) Oa A-—0. 
k=1 

(since Ly = Lr and by the definition of the stochastic integral Ito as a and o are constant ) . 


Appendix A 
Let Y, X1, X2,..., Xq be integrable square random variables defined on the same probability 
space {0, F, P} such that EX, =--- = EXg=EY=0. 
We try to find a coefficient vector b = (b,,...,6g)7 so that 
E(Y — bX, —+++ —baXa)*? = E(Y -b'X)* = min (Y aX? (Al) 
ae 


Let us denote EX =(EX),...,EXa4)™, Var(X) = [Cow X;, X;), 1,7 = 1,2,...,d] = EX XT. 


Proposition. nghieng The vector b minimizing E(Y — a? X)? is a solution of the following equation 
system : 


Var(X)b = E(XY). (A2) 
Putting U = Y —-b7X =Y —Y, with Y =b7X, then 
E(U?) = EY? —- b' E(XY) > 0. (A3) 
E(UX;) =0 for alli =1,...,d. (A4) 
EY? = EU? + EY’. (A5) 
EYY By?\'” 


p (A6) 


~ reve” (Ev) 
(p is called multiple correlation coefficient of Y relative to X). 
Proof. Suppose at first that Var(X) is a positively definite matrix. For each a € R? We have 
F(a) = E(Y — a? X)? = EY? - 2a7 E(XY) + a EX X70 (A7) 
VF(a) = -2E(XY) + 2Var(X Ja. 


OF(a) .. 
‘ = 1,2,...,d] = 2Var(X). 
It is obvious that the vector b minimizing F(a) is the unique solution of the following equation: 


VF(a) = 0 or (A2) 


N.V. Huu, VQ. Hoang / VNU Journal of Science, Mathematics - Physics 23 (2007) 143-154 153 


and in this case (A2) has the unique solution : 
b = [Var(X)] TE(XY). 
We assume now that 1 < Rank(Var(X)) =r < d. 


We denote by e), €2,...,€¢ the ortho-normal eigenvectors w.r.t. the eigenvalues j, A2,..., Ad 
of Var(X) , where A, > A2 > --- > A, > O = Ae4y = +++ = Ag and P is a orthogonal matrix with 
the columns being the eigenvectors e), €2,...,€g, then we obtain : 


Var(X) = PAP™, with A = Diag(A;, Ao, .-., Ag): 


Putting 
Z= PTX = (ek Xe X,...,e,X]', 
Z is the principle component vector of X, we have 


Var(Z) = P? Var(X)P = A = Diag(dy, Ao, ..-, Ar, 0,---, 0). 


Therefore 
EZ?,, =: = EZ3=0, so Z-4,=---=Zg=0 P-as. 
Then 
F(a) = E(Y -a?X)? = E(Y - (a? P)Z)* 
= E(Y -ajZ, —----—a5Za) 
= E(Y -afZ, ---- -—a'Z,)*. 
where 
a®? = (af,...,a3) =a! P, Var(Z,...,Z,) = Diag(Aj, A2,--+) Ar) > 0. 
According to the above result (bf,..., 0%)? minimizing E(Y — at Zi —-:- —a*Z,)? is the solution of 
mM oe 2 Hi | EZ\Y 
vee ree eee Poe | = = (A8) 
O see) Ge br] EX,Y 
or 
‘iar SF WO oe FU bt EZ,Y \ EZ\Y 
0 Ar O 0 bo | | EY) | BZY (A9) 
0 0 0 0 ea | 0 | Baa 
O ae 2 OD 2oe. 3D by 0 | \ EZaY 


with 67 ,,,...,0% arbitrary . 
Let b = (bi,..., bg)” be the solution of b7 P = b*7 , hence b = Pb* with b* being a solution of (A9). 
Then it is follows from (A9) that 


Var(Z)P7b = E(ZY) = P™E(XY) 
P™Var(X)PP™b = P? E(XY) ( since Var(Z) = P’ Var(X)P) 


or 
Var(X)b = E(XY) 
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which is (A2). Thus we have proved that (A2) has always a solution ,which solves the problem (A1). 
By (A7) , we have 


F(b) = min E(Y — re, 
= EY? — 2b7E(XY) +b? Var(X)b 
= EY* - 2b’ E(XY) + b7 E(XY) 
= EY? — bT E(XY) > 0. 
On the other hand 
EUX; = E(XiY) — E(X;b" X) =0, (A10) 


since b is a solution of (A2) and (A10) is the ith equation of the system (A2). 
It follows from (A10) that 


E(UY) =0 and EY? = E(U+Y)? = EU? + EY? + 2E(UY) = EU? + EY?. 


Remark. We can use Hilbert space method to prove the above proposition. In fact, let H be the set of 
all random variables £’s such that EE = 0, E&? < 00, then H becomes a Hilbert space with the scalar 
product (€,¢) = E&C, and with the norm ||€|| = (£€?)!/2 . Suppose that X), Xo,..., Xa, YE H, L 
is the linear manifold generated by X,, X2,...,Xq. We want to find a Y € H so that ||Y — Y|| 
minimizes, that means Y = 67 X solves the problem (A1). It is obvious that Y is defined by 


Y = Proj, Y =X andU =Y -Y € LL’. 


Therefore (Y —b7 X, X;) = 0 or E(b? X Xj) = E(X.Y) for alli =1,...,d or b? E(X7X) = E(XY) 
which is the equation (A2). The rest of the above proposition is proved similarly. 
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